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SUMMARY 


The  minimum-mass  structural  efficiency  curve  is  determined  for  sandwich- 
blade  stiffened  composite  compression  panels  subjected  to  buckling  and  strength 
constraints.  High  structural  efficiencies  are  attainable  for  this  type  of  con¬ 
struction.  A  method  of  analysis  is  presented  for  the  buckling  of  panels  of  this 
configuration  which  shows  that  buckling  of  such  panels  is  strongly  dependent  on 
the  through-the-thickness  transverse  shearing  of  the  stiffener.  Experimental 
results  are  presented  from  seven  test  specimens  made  of  layered  graphite/epoxy 
having  stiffeners  of  sandwich  construction  with  an  aluminum  honeycomb  core,  and 
there  is  good  agreement  between  theory  and  experiment. 


INTRODUCTION 

Determination  of  the  minimum  mass  of  a  stiffened  structural  panel  con¬ 
structed  of  laminated  filamentary  materials  to  carry  a  specific  compression 
load  in  the  plane  of  the  panel  is  an  important  design  problem  in  aerospace 
applications.  Stiffened  panels  having  open-section  stiffeners  are  generally 
easier  to  inspect  and  fabricate  than  panels  having  closed-section  stiffeners. 
Conventional  open-section  stiffened  panels  are  not  as  efficient,  however,  as 
closed-section  stiffened  panels  (ref.  1).  One  way  to  increase  the  performance 
of  open-section  stiffened  panels  is  by  using  sandwich-blade  stiffeners.  In 
this  configuration  the  webs  of  blade  stiffeners  are  composed  of  faces  of  high- 
strength  material  bonded  to  a  low-density  core  such  as  aluminum  honeycomb.  By 
moving  the  high-strength  material  out  from  the  blade  midsurface  the  critical 
twisting  and  local  buckling  strains  can  be  increased  without  appreciably 
increasing  the  total  structural  mass.  Some  of  the  advantage  resulting  from  mov¬ 
ing  the  high-strength  material  out  from  the  midsurface  is  offset  by  the  decrease 
in  stiffness  due  to  one  face  shearing  with  respect  to  the  other.  This  through- 
the-thickness  transverse  shearing  depends  on  the  shear  stiffness  of  the  core. 

Simple,  closed-form  column  and  local  buckling  formulas  which  treat  various 
buckling  modes  independently  are  adequate  for  the  design  of  closed-section  con¬ 
figurations,  such  as  the  hat-stiffener  configuration  (ref.  2).  More  sophisti¬ 
cated  buckling  analyses  are  required,  however,  to  predict  the  buckling  behavior 
of  open-section  configurations  due  to  their  more  complicated  stiffener  twisting 
and  rolling  behavior  (ref.  1).  For  a  conventional  stiffened  panel  acting  as  a 
column,  inplane  shearing  stiffness  of  the  webs  of  the  stiffeners  determines  the 
overall  transverse  shearing  stiffness  of  the  column.  For  a  stiffened  panel 
with  sandwich  webs  through-the-thickness  shearing  stiffness  of  the  webs  as  well 
as  inplane  shearing  stiffness  of  the  webs  must  be  considered.  In  this  regard 
finite-element  analyses  are  available  which  include  all  the  effects  needed. 

Such  analyses,  however,  typically  have  execution  times  that  are  longer  than 
desirable  for  the  many  calculations  required  in  a  structural  efficiency  study. 
Accordingly,  to  fill  this  need,  a  simplified,  efficient,  and  reasonably  accu¬ 
rate  analysis  has  been  derived  and  is  presented  in  this  paper  for  buckling  of 
laminated  sandwich-blade  stiffened  panels. 


In  addition,  the  paper  contains  a  design  procedure  based  on  the  newly 
derived  analysis  which  is  used  to  determine  the  minimum  mass  necessary  to 
carry  a  prescribed  ratio  of  compressive  load  to  panel  length  for  araphite/epcxy 
sandwich-blade  panels,  and  results  are  presented  as  structural  efficiency  plots 
and  are  compared  with  previous  results.  The  design  variables  are  the  panel 
dimensions  and  ply  thicknesses,  and  the  constraints  are  buckling  and  material 
allowable  strength.  The  validity  of  the  analytical  derivations  is  established 
by  comparison  with  a  finite-element  structural  analysis  and  with  test  results 
for  seven  typical  panels. 


L  panel  length 

Mx,My,MXy  moments  per  unit  length 
m  number  of  buckles  along  length 

N  axial  stress  resultant  of  panel 

Nx,Ny,NXy  stress  resultants 
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n 


integer 


QX,Q y  transverse  shear  stress  resultants 

tj[  thickness  of  layer  i 

u,v,w  displacements  in  x,y,z  directions,  respectively 

W  panel  mass 

x  axial  coordinate 


yj_  cross-sectional  coordinate  along  segment  i 

Bx,By  rotations  of  normal  in  web 

F  warping  constant 

Yxz'Yyz  transverse  shearing  strains 

A-j_j  Kroneker  delta,  1  if  i  =  j  or  0  if  i  ^  j 

£  axial  strain  at  buckling 

£x'ey'”Yxy  strains 

0  angle  of  twist 

Kx,Ky,KXy  curvatures 


MINIMUM-MASS  PANEL  DESIGN 
Design  Procedure 

A  computational  procedure  was  developed  for  the  design  of  minimum-mass 
sandwich-blade  stiffened  panels  loaded  in  axial  compression.  Nonlinear  mathe¬ 
matical  programing  techniques  reported  in  reference  3  are  used  in  this  design 
procedure  in  which  the  cross-sectional  dimensions  of  the  panel  are  the  design 
variables  and  buckling  and  material  strength  are  the  constraints.  The  buckling 
analysis  is  provided  by  the  present  theory  and  the  strength  analysis  by  maximum 
allowable  strain.  In  design  studies  the  panel  length  is  fixed  at  76.2  cm.  In 
the  present  analysis  it  is  assumed  that  the  compression  panel  is  simply  sup¬ 
ported  and  has  many  stiffeners  spaced  along  its  width,  so  that  the  structural 
behavior  can  be  modeled  by  a  typical  repeating  element.  The  geometry  and 
design  variables  used  to  represent  the  repeating  element  for  the  sandwich-blade 
configuration  are  shown  in  figure  1.  The  cross  section  is  composed  of  four 
characteristic  elements:  (1)  the  skin  between  two  adjacent  attached  flanges, 

(2)  the  skin  and  attached  flange  element,  (3)  the  sandwich  web,  and  (4)  the 
high-stiffness  cap.  The  sandwich  web  and  cap  constitute  the  stiffener.  Seven 
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thicknesses  and  four  width  dimensions  constitute  the  design  variables  for  the 
sandwich-blade  configuration  with  lamina  ply  orientation  angles  of  0°,  S0°, 
and  ±45°  permitted.  Orthctropic  stiffnesses  are  assumed  for  all  elements,  and 
the  additional  coupling  effects  introduced  by  nonsymmetr ic  layup  patterns  are 
ignored , 


Design  Results 

The  design  procedure  was  used  to  determine  minimum-mass  designs  required 
to  carry  a  specified  load  for  the  sandw’ich-blade  stiffened  configurations  con¬ 
structed  of  graphite/epoxy  and  aluminum-honeycomb  material.  Material  proper¬ 
ties  used  in  these  studies  are  presented  in  table  I.  Thicknesses  were  assumed 
to  be  continuous  variables  (that  is,  an  integer  number  of  plies  was  not 
required),  and  neither  thickness  nor  cross-sectional  dimension  design  variables 
were  constrained  to  minimum  limits.  To  satisfy  the  strength  requirement  a  max¬ 
imum  allowable  axial  strain  (in  direction  of  load)  of  0.008  was  specified. 

The  structural  efficiency  curve  shown  in  figure  2  was  developed  by  first 
generating  with  the  design  procedure  seven  minimum-weight  designs  in  the  lead 
index  (N/L)  range  from  0.3  MPa  tc  7  MPa.  Arbitrary  honeycomb-core  densities 
ranging  from  24  kg/nP  to  160  kg/m^  were  permitted  in  determining  the  minimum 
mass  of  a  panel  required  to  carry  a  specified  compression  load.  Transverse 
shear  modulus  properties  for  the  aluminum  honeycomb  were  taken  from  reference  4 
and  are  listed  in  table  1(b).  It  was  found  that  for  minimum  mass  the  core  den¬ 
sity  varied  from  24  kg /me  for  lightly  loaded  (N/L  =2.5  MPa)  to  160  kg/m^  for 
heavily  loaded  (N/L  =  7  MPa)  panels. 

The  structural  efficiency  curve  for  the  solid  graphite/epoxy  blade  config¬ 
uration  presented  in  figure  2  is  from  reference  1  in  which  transverse  shear 
effects  were  not  included.  The  results  presented  in  figure  2  for  the  hat- 
stiffened  configuration  are  taken  from  reference  5  and  represent  slightly 
higher  material  stiffness  properties  than  were  used  in  the  present  study. 

Comparisons  of  the  results  shown  in  figure  2  show;  the  graphite/epoxy 
sandwich-blade  configuration  to  be  approximately  30  percent  more  efficient 
than  the  graphite/epoxy  solid-blade  configuration  (that  is,  the  sandwich-blade 
configuration  requires  30  percent  less  mass  than  the  solid-blade  configuration 
to  carry  the  same  load)  and  15  percent  less  efficient  than  the  graphite/epoxy 
hat-stiffened  configuration.  Comparing  the  graphite/epoxy  sandwich-blade  con¬ 
figuration  results  with  results  presented  in  reference  6  for  aluminum-blade 
and  hat-stiffened  panels  shoves  that  the  graphite/epoxy  sandwich-blade  config¬ 
uration  provides  a  6C-percent  mass  savings  over  the  aluminum-blade  design  and 
a  50-percent  mass  savings  over  the  aluminum-hat  design. 

A  separate  investigation  was  carried  out  to  check  the  effect  on  the  struc¬ 
tural  efficiency  of  blade  web  transverse  shearing  for  the  solid  graphite/epoxy 
blade  configuration  studied  in  reference  1.  By  assuming  the  transverse  shear 
modulus  only  one-fourth  of  the  shear  modulus  (see  table  I)  and  by  assigning  0° 
oriented  graphite/epoxy  material  properties  for  the  core,  an  efficiency  curve 
was  obtained  using  the  present  approach.  The  present  results  obtained  for  the 
graphite/epoxy  solid-blade  configuration  are  only  about  3  percent  heavier  than 
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results  reported  in  reference  1  in  which  through-the-thickness  transverse  shear 
effects  were  not  included. 


BUCKLING  ANALYSES 
Description 

Accurate  finite-element  analyses  typically  have  execution  times  that  are 
longer  than  is  desirable  for  the  many  calculations  required  in  a  structural 
efficiency  study.  Consequently,  a  rapid  simplified  buckling  analysis  at  least 
an  order  of  magnitude  faster  than  available  finite-element  procedures  has  been 
developed.  This  analysis,  capable  of  addressing  the  complex  twisting  and  roll¬ 
ing  modes  of  buckling  behavior  of  sandwich-blade  stiffened  panels,  is  presented 
in  appendix  A.  The  analysis  predicts  panel  buckling  loads,  considering  a  vari¬ 
ety  of  buckling  modes  that  include  column  buckling,  stiffener  twisting  and 
rolling,  and  local  buckling,  as  well  as  the  interaction  between  these  modes. 

A  second  analysis  is  presented  in  appendix  B  and  applies  only  to  local  skin 
buckling  where  the  blade  stiffener  does  not  deform.  The  analysis  in  appendix  B 
is  a  more  accurate  representation  for  local  skin  buckling  than  that  given  in 
appendix  A. 

Buckling  analyses  developed  in  the  appendixes  consider  the  problem  of  sim¬ 
ply  supported  sandwich-blade  stiffened  wide  panels  loaded  in  axial  compression. 
It  is  assumed  in  the  development  of  these  analyses  that  the  buckling  of  a  wide 
multistiff ened  panel  may  be  analyzed  by  studying  the  buckling  response  of  a 
single  repeating  element  of  the  panel  cross  section.  The  repeating  element  is 
similar  to  that  shown  in  figure  1  except  that  the  skin  and  attached  flange  have 
a  common  center  line  as  indicated  by  the  sketch  in  appendix  A.  For  the  general 
study  in  appendix  A  the  displacement  functions  assumed  restrict  the  deformations 
to  a  mode  antisymmetric  about  the  blade  stiffener.  For  the  local  buckling  mode 
study  in  appendix  B,  however,  the  deformations  may  be  symmetric.  Results  using 
the  theory  of  appendixes  A  and  B  are  compared  later  in  this  paper  with  other 
theory  and  with  experiment. 

In  the  present  analysis  all  through-the-thickness  transverse  shearing  in 
the  web  in  the  unloaded  direction  is  completely  constrained  but  transverse 
shearing  in  the  loaded  direction  is  allowed.  A  separate  analytical  study  not 
presented  here  showTed  that  the  buckling  lead  of  the  panels  investigated  was  not 
significantly  affected  by  the  assumption  that  the  transverse  shear  stiffness  in 
the  unloaded  direction  is  infinite. 


Analytical  Comparisons 

To  determine  whether  or  not  the  present  results  are  sufficiently  accurate 
for  the  important  modes  of  buckling  and  for  various  levels  of  transverse  shear¬ 
ing  the  present  computed  approximate  results  are  compared  to  other  computed 
results  for  two  typical  panels.  These  computations  are  also  used  to  indicate 
the  significance  of  the  transverse  shearing  effects.  In  a  later  section  of 
this  paper  these  computed  results  are  compared  with  experimental  results. 
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Typical  designs  for  comparative  study.-  Two  designs,  which  are  designated 
as  moderately  and  lightly  loaded,  are  considered  which  are  near  minimum-mass 
proportions  necessary  to  satisfy  buckling  and  strength  constraints  for  a 
76 . 2-cm-long  simply  supported  panel.  The  cross-sectional  dimensions  and  lam¬ 
inate  descriptions  of  the  panel  designs  considered  are  shown  in  figure  3.  (For 
bracketed  expressions  defining  laminate  layups  in  fig.  3,  see  ref.  7.  HC  refers 
to  aluminum  honeycomb.)  The  panels  are  of  laminated  graphite/epoxy  construction 
except  for  the  web  of  the  blade  stiffener.  The  web  is  of  sandwich  construction 
with  an  aluminum  honeycomb  core  faced  with  plies  of  ±45  graphite/epcxy .  The 
high-axial-stiffness  graphite/epoxy  region  at  the  outstanding  cap  of  the  blade 
stiffener  (away  from  the  skin)  is  assumed  to  act  as  a  beam  attached  to  the  blade 
web.  Although  the  panel  configuration  shown  in  figure  3  has  eccentricities  in 
the  skin  to  permit  a  smooth  outer  surface,  these  eccentricities  of  the  attached 
flange  are  ignored  in  the  analysis.  The  panel  designated  as  moderately  loaded 
is  designed  to  carry  1.36  MN/m  when  the  web  is  constructed  with  a  honeycomb  core 
of  density  160  kg/m3  and  to  carry  1.15  MN/m  when  constructed  with  a  48  kg/m3 
core  density.  The  lightly  loaded  panel  is  designed  to  carry  0.53  MN/m  when  con¬ 
structed  with  a  48  kg/m-5  honeycomb-core  density.  The  buckling  strains,  which 
are  based  on  the  derivations  given  in  appendixes  A  and  B  for  the  two  different 
designs,  and  from  which  comparisons  are  made  between  the  different  modes  of 
buckling,  are  presented  in  figures  4  and  5  as  a  function  of  buckle  length,  L/m. 

Results  for  moderately  and  lightly  loaded  designs.-  Results  for  the  moder¬ 
ately  loaded  design  plotted  in  figure  4  show,  for  the  interim  region  of  buckle 
length  (from  7  cm  to  70  cm),  that  the  buckling  strain  is  strongly  dependent  on 
core  density  and,  therefore,  strongly  dependent  on  the  through-the-th ickness 
transverse  shear  stiffness  of  the  stiffener:  for  example,  for  a  30-cm-long 
buckle  length  the  48  kg/m3  density  core  reduces  the  buckling  strain  by  40  per¬ 
cent  from  the  solution  for  the  infinite  transverse  shear  stiffness  case.  Simi¬ 
larly,  for  the  same  buckle  length  a  panel  with  a  160  kg/m3  core  reduces  the 
buckling  strain  by  18  percent.  In  this  interim  range  buckling  deformations 
involve  twisting  and  rolling  of  the  stiffener.  In  the  range  beyond  70  cm  the 
buckling  strain  is  not  as  strongly  dependent  on  cere  density.  In  the  range  of 
longer  buckle  lengths,  column  effects  couple  with  some  twisting  of  the  stiff¬ 
ener;  then,  with  still  longer  lengths  column  buckling  couples  with  conventional 
column  transverse  shearing  (shearing  through  the  depth  of  the  stiffener  but  not 
through  the  stiffener  thickness)  .  In  the  range  below  7  cm  the  buckling  mode  is 
local  with  only  the  skin  and  the  attached  flange  deforming,  and  the  present 
results  are  calculated  using  the  analysis  of  appendix  B.  Of  course,  with  the 
stiffener  not  deforming  there  is  no  effect  of  core  density  of  the  stiffener  in 
this  range. 

Results  for  the  lightly  loaded  design  plotted  in  figure  5  have  ranges  sim¬ 
ilar  to  those  indicated  for  the  moderately  loaded  design.  Again  for  the  interim 
region  of  buckle  length  (from  7  cm  to  70  cm)  the  buckling  strain  is  strongly 
dependent  on  core  density  and,  therefore,  on  through-the-thickness  transverse 
shearing  of  the  stiffener:  for  example,  for  a  30-cm-long  panel  the  48  kg/m3 
density  core  reduces  the  buckling  strain  by  approximately  25  percent  from  the 
solution  for  the  infinite  transverse  shear  stiffness  case. 

In  each  of  the  design  studies  the  honeycomb  core  was  oriented  with  the 
stiffer  direction  parallel  to  the  lead  axis  of  the  panel  (that  is,  along  the 
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length  of  the  panel;  see  table  I) .  Reversing  this  orientation  was  found  to 
decrease  the  calculated  buckling  strain  in  the  interim  stiffener-twisting¬ 
mode/stiffener-rolling-mode  range  by  approximately  15  percent.  To  achieve 
maximum  load,  therefore,  it  is  recommended  that  the  core  be  oriented  so  that 
its  stiffest  direction  is  along  the  load  axis. 

Comparisons  with  BUCLASP2.-  When  the  honeycomb  core  is  very  dense,  trans¬ 
verse  shearing  is  suppressed  and  the  blade  stiffener  web  may  be  treated  as  a 
conventional  plate.  An  exact-plate-analysis  BUCLASP2  is  available  in  refer¬ 
ence  8  for  the  buckling  of  an  assembly  of  conventional  plates  and  beams  which 
includes  inplane  shearing  but  does  not  allow  through- the- thickness  transverse 
shearing  deformations.  Comparison  of  present  results  and  the  exact  results 
using  the  method  of  reference  8  (figs.  4  and  5)  shows  that  the  present  results 
for  infinite  core  density  follow  the  same  trends  as  exact  results  and  agree  to 
within  a  few  percent  throughout  most  of  the  buckle  length  range.  In  the  range 
below  7  cm  when  local  buckling  occurs  there  is  better  agreement  for  moderately 
loaded  panels  (fig.  4)  than  for  lightly  loaded  panels  (fig.  5)  because  the 
inplane  deformation  of  the  stiffener  that  occurs  when  the  faces  are  very  thin 
is  not  modeled  by  the  present  local  buckling  analysis. 

Comparisons  with  NASTRAN-*-.-  In  order  to  assess  the  accuracy  of  the  present 
analysis  to  account  for  through- the- thickness  transverse  shear  effects  compari¬ 
sons  have  been  made  using  the  finite-element  analysis  code  NASTRAN  (ref.  9) , 

Use  of  standard  QUAD1  elements  in  NASTRAN  allowed  for  orthotropic  plate  proper¬ 
ties  and  isotropic  transverse  shearing.  Comparison  of  results  in  figure  4  shows 
that  the  methods  of  references  8  and  9  agree  for  infinite  core  density.  For 
finite  core  density  there  is  also  very  good  agreement  between  present  results 
and  the  results  obtained  using  NASTRAN. 


EXPERIMENTAL  RESULTS 
Specimen  Description  and  Test  Technique 

The  moderately  and  lightly  loaded  sandwich-blade  stiffened  compression 
panel  designs  studied  analytically  in  the  previous  section  were  also  studied 
experimentally.  A  summary  of  design  load  and  mass  of  the  panels  is  given  in 
table  II,  and  the  cross-sectional  dimensions  and  laminate  descriptions  for 
these  two  designs  are  presented  in  figure  3.  The  specimens  were  tested  in  a 
two-stiffener-wide  configuration  as  shown  in  figure  3  with  side  edges  unsup¬ 
ported  during  the  test.  The  specimen  ends  were  potted  in  epoxy  approximately 
2.5  cm  thick,  were  ground  flat  and  parallel,  and  were  tested  "flat  ended"  for 
a  nearly  clamped-end  test  condition.  Panels  were  tested  in  lengths  ranging 
from  17.9  cm  to  76.2  cm. 

Specimens  were  constructed  using  7.6-cm-wide  preimpregnated  tape  of 
graphite  in  epoxy  resin  and  aluminum  honeycomb.  Mechanical  properties  for 
these  materials  are  presented  in  table  I.  Moderately  loaded  specimens  were 


^NASTRAN:  Registered  trademark  of  the  National  Aeronautics  and  Space 
Administration. 
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constructed  using  a  honeycomb  density  of  either  48  kg/m3  or  160  kg/m3.  Lightly 
loaded  specimens  used  a  48  kg/m 3  density  core.  The  concentration  of  0°  ori¬ 
ented  plies  at  the  stiffener  cap  was  laid  up  but  not  cured  in  a  separate  opera¬ 
tion  and  was  transferred  to  the  stiffener  as  a  unit.  Epoxy  was  used  to  fill  the 
honeycomb  core  for  approximately  0.38  cm  at  the  junctures  where  the  honeycomb 
meets  the  skin  and  cap.  A  separate  adhesive  layer  was  positioned  in  the  web 
between  the  honeycomb  and  the  graphite  epoxy.  The  entire  assembly  was  vacuum- 
bagged  and  cured  in  an  autoclave.  Adhesive  and  epoxy  materials  for  the  core 
were  not  accounted  for  in  the  calculations  and  resulted  in  panels  that  were 
heavier  than  was  estimated.  Representative  measurements,  which  are  presented 
in  table  II,  indicated  that  the  panels  were  5  to  33  percent  heavier  than  calcu¬ 
lated.  The  lightly  loaded  panels  had  more  adhesive  materials  for  the  core  than 
the  moderately  loaded  panels  on  the  basis  of  percent  of  total  weight.  Other 
measurements  of  the  panels  indicated  very  close  to  nominal  dimensions.  Experi¬ 
mental  buckling  strains  were  obtained  from  strain  reversal  data,  and  the  mode 
shapes  were  determined  from  displacement-gage  data  and  moire  fringe  patterns 
(ref.  10) . 


Experimental-Analytical  Comparison 

A  comparison  of  experimental  buckling  strains  and  mode  shapes  between  the 
results  obtained  using  NASTRAN  and  the  present  theory  is  presented  in  table  III. 
Graphical  comparisons  of  the  buckling  strain  as  a  function  of  the  panel  length 
are  presented  in  figures  6  and  7.  Analytical  results  using  the  analysis  devel¬ 
oped  herein  are  for  a  simply  supported  infinitely  wide  multistiffener  configura¬ 
tion  while  the  NASTRAN  results  and  experimental  data  are  for  the  clamped  two- 
stiffener  configurations  shown  in  figure  1. 

Moderately  loaded  design.-  Good  agreement  for  buckling  strains  was  achieved 
between  both  theories  and  experiments  as  shown  in  figure  6  and  table  III.  Each 
of  the  48  kg/m3  honeycomb  core  density  specimens  buckled  in  a  mode  involving 
stiffener  twisting  and  buckling  of  the  skin  as  predicted  by  theory.  Increasing 
the  honeycomb  core  density  from  48  kg/m3  to  160  kg/m3  increased  the  buckling 
strain  by  approximately  40  percent  and  changed  the  mode  shape  from  m  =  2  to 
m  =  1  as  predicted  by  theory.  The  buckling  deformation  pattern  for  the  skin 
segments  of  the  test  panels  may  be  interpreted  using  the  moire  fringe  photo¬ 
graphs  shown  in  figure  8.  Back-to-back  strain-gage  data  and  deflection-gage 
data,  the  results  of  which  are  not  presented,  helped  to  indicate  that  stiffener 
twisting  and  rolling  were  occurring.  In  recording  mode  shapes  the  buckle  of  the 
skin  between  stiffeners  at  the  panel  ends  has  been  ignored  for  long  wavelength 
buckling  since  this  condition  is  peculiar  to  the  clamped-end  test  restraint 
(ref.  1)  and  the  corresponding  clamped-end  conditions  imposed  on  the  NASTRAN 
solution.  Good  correlation  for  mode  shape  was  obtained  between  NASTRAN  and 
experimental  results. 

Lightly  loaded  design.-  Good  agreement  for  buckling  strains  was  obtained 
between  the  results  from  the  present  theory  and  the  experimental  results 
obtained  for  both  of  the  40-cm-long  lightly  loaded  test  panels  as  shown  in 
figure  7  and  table  III.  However,  both  panels  buckled  in  a  local  skin  mode 
(m  =  7)  rather  than  in  a  stiffener-twisting/skin-buckling  mode  (m  =  2)  as  pre¬ 
dicted  by  the  present  theory.  (See  table  III.)  A  moire  fringe  photograph  of 
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the  buckle  pattern  is  presented  below  the  curve  in  figure  7.  It  is  not  clear 
why  the  experimental  mode  shapes  do  not  agree  with  theory  for  the  lightly  loaded 
design  while  they  did  agree  with  theory  for  the  moderately  loaded  design.  From 
a  study  of  figure  4,  however,  the  minimum  strain  for  short  wavelength  buckling 
(L/m  <  7  cm)  is  close  to  the  minimum  strain  for  interim  wavelength  buckling 
(7  cm  <  L/m  <  70  cm)  for  the  lightly  loaded  panel.  Thus,  the  determination  of 
the  mode  shape  is  sensitive  to  the  detailed  test  conditions. 

The  ends  of  the  specimens  were  all  ground  flat  and  parallel  and  tested 
"flat-ended"  so  that  they  are  essentially  clamped.  The  present  theory  assumes 
the  ends  are  simply  supported.  In  the  interim  range  of  buckle  length  where 
twisting  occurs,  however,  the  results  using  the  present  theory  agree  well  with 
clamped  experimental  and  NASTRAN  results  (fig.  4) .  Because  of  the  good  agree¬ 
ment  between  clamped  test  results  and  both  simple-support  and  clamped  theory, 
it  is  concluded  that  detailed  boundary  conditions  beyond  requiring  zero  deflec¬ 
tions  are  apparently  not  important  for  the  buckling  of  sandwich-blade  stiffened 
panels  where  twisting  occurs  and  where  through-the-thickness  transverse  shearing 
is  important. 


CONCLUDING  REMARKS 

A  minimum-mass  study  is  presented  for  sandwich-blade  stiffened  compression 
panels  subject  to  buckling  and  strength  constraints.  This  study  indicates  that 
high  structural  efficiencies  are  attainable  for  open-section  stiffened  panels. 
The  graphite/ epoxy  sandwich-blade  stiffened  configuration  is  approximately 
30  percent  more  efficient  than  a  graphite/epoxy  nonsandwich  configuration, 

50  percent  more  efficient  than  an  aluminum-hat  stiffened  configuration,  and 
60  percent  more  efficient  than  an  aluminum-blade  stiffened  configuration. 
Analysis  suitable  for  use  in  automated  design  is  presented  for  the  buckling  of 
sandwich-blade  stiffened  compression  panels.  The  buckling  strain  of  sandwich- 
blade  stiffened  panels  depends  strongly  on  the  transverse  shear  stiffness  of 
the  blade  web  core.  The  results  of  the  buckling  analysis  which  includes  trans¬ 
verse  shear  deformations  agree  well  with  finite-element  results.  Buckling 
experiments  on  clamped,  one-bay,  graphite/epoxy  panels  having  two  sandwich- 
blade  stiffeners  with  aluminum  honeycomb  core  are  also  presented,  and  the 
results  of  the  buckling  analysis  also  agree  with  these  experiments.  Conse¬ 
quently,  the  analysis  and  especially  the  inclusion  in  the  analysis  of  the 
effects  of  transverse  shearing  of  the  blade  web  core  and  the  assumption  that 
the  panel  consists  of  an  infinite  number  of  repeating  bays  are  validated. 

The  theoretical  results  using  simple-support  boundary  conditions  agree  with 
clamped  theoretical  and  experimental  results  when  twisting  occurs  and  through- 
the-thickness  transverse  shearing  is  important. 


Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Hampton,  VA  23665 
August  3,  1978 
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APPENDIX  A 


BUCKLING  ANALYSIS  OF  SANDWICH-BLADE  STIFFENED  COMPRESSION  PANELS 

Buckling  of  wide,  simply  supported,  sandwich-blade  stiffened  panels 
loaded  in  axial  compression  is  analyzed.  The  analysis  predicts  the  buckling 
load  considering  a  variety  of  buckling  inodes  including  column  buckling,  stiff¬ 
ener  twisting  and  rolling,  local  buckling,  and  interaction  between  the  modes. 
It  is  assumed  that  the  buckling  of  a  wide  multistif fened  panel  may  be  analyzed 
by  studying  the  repeating  element  of  the  panel  cross  section,  which  is  shown 
in  the  following  sketch: 


The  model  is  composed  of 
flange  plate  segments  (?) 
and  a  sandwich  web  plate 


six  basic  segments:  two  skin  (?)  and  two  attached 
located  symmetrically  and  having  a  common  center  line, 
segment  (?)  connected  to  a  beam  (?)  forming  a  blade 


stiffener.  The  skin,  flange,  web,  and  beam  segments  are  numbered  1  to  4, 
respectively,  and  the  coordinate  system  is  indicated  in  the  sketch.  The  small 
error  associated  with  analyzing  a  wide  panel  by  studying  the  response  of  a  sin¬ 
gle  repeating  element  of  the  panel  is  discussed  in  reference  1. 


Expressions  are  obtained  for  the  buckling  strain  e  in  this  appendix  by 
the  method  of  virtual  work  in  conjunction  with  assumed  displacements.  The  vir¬ 
tual  work  of  extension,  bending,  and  the  applied  strain  e  at  buckling  for  the 
various  plate  segments  may  be  expressed  as 
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pfc>i  p  L 

6tt  —  |  I  Ny6£j^  +  NyS£y  +  NXy6yXy  ^■‘x'^X  ^y'-’^y 

<vn 


"0  '-'O 


+  Mxy6Kxy  -  el  An  -  —  j(w>x6w/X  +  v(X6v(X)  dx  dy  (Al) 


where  the  N's  and  M's  are  the  stress  resultants  and  moments;  the  e 's,  YXy' 
and  k's  are  the  strains  and  curvatures;  and  w  and  v  are  deformations  that 
occur  during  buckling.  The  axial  length  of  the  panel  is  L,  and  the  width  of 
each  segment  is  b-j_.  The  strains  and  curvatures  are  given  in  terms  of  the  dis¬ 
placements  as  follows: 


ex  ~  u,x 


ey  =  v»y 


Yxy  =  u,  y  +  vfX 


Kx  ~  w,xx 


Ky  -  W,  yy 


^xy  _  2wfXy 


The  stress-strain  and  moment-curvature  relations  are: 


Nx  =  Allex  +  A12ey 
Ny  =  A22ey  +  A12ex 
Nxy  =  A66Yxy 


Mx  =  DllKx  +  D12Ky 
My  =  D22Ky  +  D12kx 

MXy  =  Dge^xy 


In  order  to  account  for  transverse  shearing  of  the  blade  stiffener,  the 
following  work  of  the  shearing  forces  must  be  added  to  the  virtual  work  expres¬ 
sion  for  the  blade  web  segment 


b3  rL 


r°  3  r 

'-'0  ^0 


(QxfiYxz  +  Qy^Yyz)  dx  dy 


For  the  web  segment  the  curvatures  are  redefined  to  be 


Kx  -  ^x,x 


K„  =  B..  .. 
y  y  i  y 


Kxy  ~  ^x,y  +  ^y,x 
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The  additional  strain-displacement  and  stress-strain  relations  for  the  web  are 


Yxz  _  ^x  w,x 

Cx  = 

\ 

(A6) 

YyZ  =  By  -  Wfy 

Qy  -  ^55Yyz  j 

For  the  beam  the  virtual  work  is 

EIlw,xx^w,xx  +  ei2v,xx^v,xx  +  E^®,xx^®,xx  +  GJ®,x^®,x  +  EAu,x<^u,x 


-  £EA 


x 


+  v  V<5v  v 

/A  f  X 


dx 


(A7) 


w'here  EIj_  is  the  beam  bending  stiffness  about  the  >’3  axis,  El 2  the  beam 
bending  stiffness  about  an  axis  in  the  cross-section  perpendicular  to  the  y 3 
axis,  ET  is  the  warping  stiffness,  GJ  the  torsional  stiffness,  EA  the 
extensional  stiffness,  Ip  the  polar  moment  of  inertia,  and  A  the  cross- 
sectional  area  of  the  beam. 

The  total  virtual  work  is  equal  to  the  sum  of  the  virtual  work  for  each  of 
the  six  segments.  By  substituting  the  expressions  (A2),  (A3),  (A5),  and  (A6) 

for  stress  resultants  and  moments  and  the  strains  and  the  curvatures,  the  vir¬ 
tual  work  equation  can  be  expressed  in  terms  of  the  displacements  u,  v,  and 
w  for  plate  segments  in  the  skin  and  in  terms  of  u,  v,  w,  Bx,  and  By  for 
the  web  plate  segment  in  the  blade  stiffener.  By  requiring  that  the  beam  be 
attached  firmly  to  the  end  of  the  web  of  the  blade  stiffener,  the  u,  v,  w, 
and  6  deformations  of  the  beam  are  given  in  terms  of  the  deformations  of  the 
web.  A  characteristic  value  problem  for  the  lowest  applied  strain  e  (the 
buckling  strain)  is  obtained  by  considering  particular  expressions  for  the  dis¬ 
placements  that  satisfy  geometric  boundary  and  continuity  conditions.  The  dis¬ 
placements  for  each  of  the  plate  segments  that  have  been  chosen  have  the  follow¬ 
ing  form  with  the  following  arbitrary  displacement  coefficients  U,  U3,  V,  W, 
W3,  w30'  x3'  x30'  an<^  Y'  plus  arbitrary  parameters  A  and  y: 

Segment  Displacements 


u  =  U  cos  - 

L 

1  v  =  0  ^  (A8) 

/  2b!  TTyp 

w  =  -IV  +  boW  +  -  W  cos  - 

\  *  2bl 
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Segment 


Displacements 


u  =  U  cos 


v  =  0 


w  =  - (V  +  Wb2  -  Wy2)  sin 


u  =  (U  +  U3Y3)  cos 


v  =  V  sin 


XY3  f  YY3\  mTTx 

w  =  W3  sin  -  +  W30II  -  cos  -  1  sin  - 

b3  \  b3  /  L 


(A10) 


/  XY3  YY3\  mTTx 

3X  =  X  +  X3  sin  -  +  X30  cos  -  cos  - 

\  b3  b3  /  L 


By  =  W  COS 


Xy3  yy3 

-  +  Y  sin  - 


The  displacements  of  the  beam  are  given  in  terms  of  the  displacements  of 
segment  3  according  to 

u  =  u(b3)  -  (b4  -  b3)vf  x(b3)"] 


v  =  v(b3) 

w  =  w(b3)  +  (b4  -  b3)3y(b3) 


(All) 


These  displacements  are  chcsen  so  that  they  satisfy  geometric  continuity  between 
segments,  simple-support  boundary  conditions  at  x  =  0,L,  and  have  zero  slope 
at  Yi  =  0  and  the  corresponding  point  on  the  other  side  of  the  repeating  ele¬ 
ment.  Displacement  functions  are  chosen  in  the  y  direction  that  give  reason¬ 
able  results  for  buckling  in  the  local,  stiffener-twisting,  and  column  modes. 
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By  assuming  the  expressions  given  for  w  in  segments  1  and  2,  accuracy  in  the 
local  (skin-only)  mode  has  been  sacrificed  for  simplicity,  A  more  accurate 
solution  for  local  skin  buckling  is  presented  in  appendix  B. 

Many  of  the  displacements  assumed  for  the  general  buckling  analysis  pre¬ 
sented  here  are  zero,  constant,  or  linear  in  the  y  direction.  The  skin- 
deflection  expression  has  a  trigonometric  function  in  the  y  direction  to 
allow  local  buckling  and  twisting.  The  blade  displacement  expressions  have  two 
trigonometric  functions  in  the  y  direction  plus  arbitrary  parameters  A  and 
Y  to  allow  for  rolling  as  well  as  local  buckling  and  twisting. 

Because  of  the  way  the  panels  are  constructed,  it  is  reasonable  to  expect 
that  in  segment  3 

By  -  wfy  =  0  (A12) 

at  y3  =  0,b3.  As  a  consequence  of  this  assumption 


(A13) 


and,  therefore,  in  effect  it  is  assumed  that  the  transverse  shearing  strain 
YyZ  is  zero  in  the  web.  The  solution  is  then  independent  of  A55,  and  the 
number  of  unknown  displacement  coefficients  is  reduced  from  10  to  8. 

The  stability  equations  that  result  from  the  method  of  virtual  work  for 
this  problem  can  be  identified  from 


a12 

a22 


Symmetric 


a13 

a23 

a33  “  eb33 


0 

0 

a34  “  eb34 
844  -  eb44 


0 

0 

0 


0  0  0 
0  0  0 
0  0  0 


a45  _  eb45 

a46 

a47 

a48 

w3 

a55  "  eb55 

a56 

a57 

a58 

w30 

a66 

a67 

a68 

:  X 

a77 

a78 

*30 

ru 

CO 

CO 

*3 

0 


(A14) 
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wherein  what  follows  A  =  A^q  -  {a12/a22)  and  subscripts  have  been  added  to  A, 
A^j,  and  D^j  to  indicate  the  segment  represented: 


nm  \ 

all  =  l~j  (2Alllbl  +  2A112b2  +  A113b3  +  EA) 


(m\2f1  o  \ 

a12  =  (“  A113b32  +  mh3j 


m\ 

a13  =  (h4  ”  b3)EA 


a22  =  f  (-  A113b33  +  EAb3'2)  +  A663b3 


mu  /  nm  \ 

a23  =  ~  A663b3  "  (  ~  I  (EAh>3 )  (b4  -  b3) 


l33 


=  2{~]  [Dlllbl  +  D112b2  +  ei2  +  EA(b4  -  b3)  2J  +  f— j  A663b3 


/  nor 

b33  =  (  —  )  <A3b3  +  2Albl  +  2A2b2  +  EA) 


X  (m  \ 

034  ■  2d(rj 


blb2  + 


(?) 


+  T  D112b22  \  +  It)  d121  I 


X  fm 

b34  =  2  ' fii1 


'  2bi  \2 


blb2  +  \~) 


1  -  9 

+  -  A2b22 


nm 


a44  “ 


1  (  cos  X  sin  X\ 

—  A44b3  1 1  -  - — - j  +  EI]_  sin2  X 


+  l7,Fl 


X  /nm 

+  2  —  — I  (b4  -  b3)  sin  X  cos  X 
b3\L 
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+  2  —  EA  (b4 
b3 


1  - 


cos 


sin 


-  b3)  sin  X  cos  X 


EA  sin^  X 


/  m  tt  \  ^  /mTr\^  X  y 

a45  =  —  I  a44f3  +  —  El j  sin  X(1  -  cos  y)  +  —  F4  +  — -  F5 

\L  /  \L  /  b3  b3 


b45  = 


X 

+  EA  sin  X(1  -  cos  y)  +  —  EA  (b4  -  b3) (1  -  cos  y)  cos  X 

b3 


F6 


mTT  b3 

a46  =  —  — (!  ~  cos  X) 
L  A 


d123  “  A4 4 


a47 


L 


mTT  X 

+  7“  F8  ~  D223 

L  bq 


X  Y 

—  —  d663fI0 
b3  b3 


a55  - 


2 

a44f11 


+ 


4 

I  EI^  (1  -  cos  Y) 


2 


,  „  ,4  ,  ,2 

y  mv  \  y  \ 

+  2  —  El j  —  (b4  -  b3 )  sin  Y  (1  -  cos  Y)  +  —  F12 

b3  \L  /  \b3 / 


>55 


v 

a3f11  +  EA(1  -  cos  Y)2  +2  —  EA (b4  -  b3)  sin  Y (1  “  CO£  Y) 

b3 


/  \2 
y  \ 

- —  EA 

2  IP 

(b4  -  b3)  +  — 

,21 
sin^  Y 

\b3/ 

A 

1 
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a56 


a57 


a58  = 


a66  = 


a67  = 


a68  = 


a77  = 


a78  = 


a88  " 


IWT 

L 


A44  b3 


siny\  Y 

-  +  d123  —  sin  Y 

Y  /  b3 


nm 

-  —  A44F3 

Lj 


iwt 

— -  d123f9 

Lj 


\y  \nfrr 


IWT 

-  —  A44F13 

Li 


Y 


j  D113b3  +  A44b3 


/  irtir  b3 

(  —J  d113  +  A44  “(1  ~  cos 


X) 


+  A44 


b3 

—  sin  Y 
Y 


+  A44 


F7 


+ 


D113  +  a44 


Fq 


a44 


f15 


d663fI0 
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'mif  j  „  /2blW  bl\  b2 

~  \  Dm  bib22  +  2{~r  ! b2  +  —  I  +  Dii2  — r 


/m7r\^  (  bl 

+  4  i  7 —  J  D121  b2  +  ~ 


/  7T  \  /m  TT  \  f  bl  \  /  m  7T  \  /X  \ 

+  77~  D221bl  +  8  ~  !d661  T"  +  D662b2  j  +  !  ~J  D663F8  +  7~  D223F7 

y 2bi  /  \L  j  y  2  j  \L  /  \^3/ 


r  i^\  -  /^v3  /m5T\2 

+  Nr  ,b<  - b31'  *  t)  +GJ  r 


ccs2  X 


/mirr  I  _  .  /2blW  bl\  -  b2' 

f2  =  2  —  \  A1  blfc22  +  2  - j  b2  +  ■  J  +  a2  - 

\L  j  ]  \  TT  /  \  4  /  3 


mTi\  0  ,  P 

+  EA( —  j  ccs2  X  (fc>4  -  b3)2  +  - — 


b3  1  -  cos  (X  -  y)  1  -  cos  (X  +  y)  b3 
F-j  =  - - +  -  +  — (1  -  cos  X) 

2  X  -  y  X  +  y  X 


EA^—  I  (b4  -  b3 )  (1  -  cos  Y)  cos  X 


/mTT  \  X 

EI][ — j  (b4  -  b2)  sin  Y  sin  X  +  —  Fj-y 


/itflT  V  X 

EA  —  (b4  -  b'3 )  sin  y  sin  X  +  —  Fj_g 

\L  /  b3 


b3 /  cos  X  sin  X 

r1- — i — 


b3  /  cos  X  sin  X 


F»  -  !  1  + 
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Fg  = 


b3 


1  -  cos  (X  -  y)  1  -  cos  (X  +  y) 

—  +  — 


X  -  Y 


X  +  y 


'10 


1  -  cos  (X  -  Y)  1  -  cos  (X  +  y) 
—  +  — 


X  -  Y 


b3  /  cos  Y  sin  Y 

m  -  r  i  ♦  - 


X  +  y 


b3 


,  +  b3  -  2  —  sin  Y 

y  y 


iitTT  \2  /  Y  \2 

■12  =  l  “  )  D663F16  +  f  —  j  D223F15 


EIi(r) ,b«  - b3)2  +  Er(r)  +  GJ(r) 


b3  /  cos  Y  sin  Y  \  b3 
f13  = - 1  +  -  +  —  sin  Y 

2  V  y  1  y 


mTT  Y 

r14  =  ~  “  ~~  (D123F15  +  d663f16) 

L  b3 


°3  cos  Y  sin  Y 
f15  =  T" U  +  ‘  ~ - 

*  \  y  1 


sin2  Y 


D3  /  cos  Y  sin  Y 

Fi6  =  ~  i - : - 


mil  \  Xy 

,  r  ,  d663f10 - 7  d223f9 

\L  /  b32 


F17  =|  — 


4- 

1 

|  ElJ 

( mTT  ^ 

i 

l  ' 

_ 

vL  ) 

4 

nrrr 

+  GJf 


sin  Y  cos  X 
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F18 


sin  Y  cos  X 


The  stability  equations  (A14)  are  a  homogeneous  set.  To  obtain  a  nontriv 
ial  solution  for  the  buckling  strain  the  determinant  of  the  matrix  must  vanish 
The  buckling  strain  is  the  lowest  value  of  the  applied  strain  e  that  satis¬ 
fies  this  equation  for  various  values  of  X,  y,  and  m.  As  in  reference  1, 
suggested  trial  values  are  X  =  O.Itt,  0.3tt,  0.5tt,  0.7tt,  and  y  =  X  +  0.5tt. 
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LOCAL  BUCKLING  OF  SANDWICH-BLADE  STIFFENED  COMPRESSION  PANELS 

An  analysis  is  developed  for  the  skin-only  local  buckling  of  wide,  simply 
supported,  sandwich-blade  stiffened  panels  loaded  in  axial  compression.  The 
method  of  analysis  used  is  the  method  of  virtual  work  in  conjunction  with  cen¬ 
tral  finite  differences.  For  this  mode  of  buckling  the  blade  stiffener  does 
not  deform  but  serves  to  clamp  the  skin  at  the  point  of  attachment.  A  sketch 
of  the  model  studied  is  shown  here: 


|  « - hi - -H  < - b2 

_ i - 

- y 


For  the  mode  considered,  the  slope  at  y  =  0  is  zero.  For  this  model  stretch¬ 
ing  and  bending  are  uncoupled  and,  for  buckling,  only  bending  must  be  consid¬ 
ered.  The  virtual  work  during  buckling  is 


cbi+b2  rL 

'  & 

= 

Mx6kx  +  MyfiKy  +  Mxy6Kxy  -  el 

AH  "  7 

!w  v6  w  v 

f  X  f  X 

O 

o 

_ 

y  A22  / 

dx  dy 


where 


(Bl) 


=  -w 


XX 


~N 

Mx  =  Dll^x  +  D12K'y 


Ky  ~  w,  yy 
Kxy  =  ~2wfXy 


My  =  D22<y  +  D12kx  \ 
Mxy  =  D66icxy 


(B2) 


nrrrx 

Letting  w  =  W (y)  sin  -  in  eauations  (Bl)  and  (B2)  and  integrating  over  x 

L 

yields 


d22w" 


D 


12 


It 


(B3) 
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where  A  =  A^q  -  (aj^/a.^)*  In  equation  (B3)  the  derivatives  are  approximated 
using  central  differences  as  indicated  in  the  following,  and  the  trapezoidal 
rule  is  used  for  integration: 

wn+l  _  2wn  +  wn-l 

(W")n  =  -  (B4) 

h2 


wn+l  “  wn 

(W)  2  =  -  (B5 ) 

n4—  h 

2 


where  h  is  the  uniform  mesh  spacing. 


Equation  (B3) 


becomes 


bn(6wn+1 
an<Swn  +  - 


26wn  +  6wn_2 ) 


h2 


An0  +  AnN 


N-l 


n=0 


where,  away  from  the  boundaries: 


(B6 ) 


a 


n 


wn+l 
D12n  - 


2Wn  +  Wn_2 
h2 


b 


n 


wn+l 
D22n  ~ 


2Wn  +  Wn_2 
h2 


(B7 ) 


c  1 


n 


2 


wn+l  “  wn 


h 


J 


and  is  the  Kroneker  delta  which  equals  zero  if  the  subscripts  are  differ¬ 

ent  and”one  if  the  subscripts  are  the  same,  and  D£j  and  A  have  the  sub- 
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script  n  that  refers  to  nodal  points.  The  expression  (B6)  for  the  virtual 
work  depends  on  the  displacements  Wn  specified  at  the  nodal  points  within  and 
on  the  boundaries  from  n  =  0  to  N  and  at  the  fictitious  points  -1  and  N+l 
outside  the  boundaries.  To  apply  the  method  of  virtual  work  it  is  necessary 
that  the  displacements  Wn  and  the  virtual  displacements  6wn  satisfy  geo¬ 
metric  boundary  conditions.  The  boundary  conditions  specified  are 

Zero  slope  at  y  =  0  or,  therefore,  W_p  =  Wp 

Zero  deflection  at  y  =  bp  +  b2  or,  therefore,  =  0 

Zero  slope  at  y  =  bp  +  b2  or,  therefore,  WN+p  =  V%_p 

These  conditions  affect  the  a's,  b's,  and  c's,  so  that  the  following  defini¬ 
tions  must  be  used  instead  of  the  corresponding  general  expressions: 


a0 


Wi  -  Wq 

2d120  - - - £A0W0 

h2 


A 


Wp  —  Wq 

b0  =  2d220  - - - 

h2 


aN-l  “ 


ran 


'mrr 


\  2 


/-2WN_p  +  WN_2\ 


(dH)n-i1  ~  J  WN-1  “  <d12)n-i1 


+  £AN_pWN_p 


>  (B8 ) 


bN-l  =  <d22)n-1  “ 


2wN-1  +  %-2 


(m 

~  (D12)N-ll  “  j  wN-l 


<d66)  1 
N — 
2 


%-l 

h 


%-l 

hN  =  2d22N  - ~" 

h2 
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The  virtual  work  become; 


6tt  = 


2  !  \  2 


Lh!  (  ao  bl  ~  b0  cl/2 


6w0  +  ap  + 


b2  -  2bp  +  b0  C3/2  -  ci/ 


cil'ci3 
bN  ~  2bN-l  +  bN-2  ^  2  N  2 

+  ^  ♦  - - ; - -/6K»- 1 


Setting  the  virtual  work  equal  to  zero  and,  therefore,  requiring  that  the 
coefficients  of  the  virtual  displacements  vanish  leads  to  the  stability  equa¬ 
tions  (for  example,  with  N  =  8) : 


all  a12 

a13 

0 

0 

0 

0 

0 

w0 

|  a22 

a23 

a24 

0 

0 

0 

0 

w2 

1 

j 

a33 

a34 

a35 

0 

0 

0 

w3 

i 

j 

a44 

a45 

a46 

0 

0 

w4 

a55 

e56 

a57 

0 

w5 

a66 

a67 

a68 

k 

Symmetric 

a77 

a78 

w7 

Q) 

CO 

CO 

1 _ 

!w8 

(BIO) 


where 


1  f  urn 
all  =  7  D110 — 


r)  +  (“121  *  0120  ♦  «66ij(“)  b  *  <20221  +  B220)L  -  £  K0^ 
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a22  =  Dlll 


+  4  Dm 


1  _  ( m7r\2 

+  (D222  +  4d221  +  2d220)‘~  “  eAl  — j 

h4  \L  / 


a33  =  d112 


d122 


+  (D223  + 


1 

4D222  +  D22l) - £A2 

h4 


a77  =  d116 


+ 


/ 


D126 


D  13 
66 — 
2 


1 

h2 


1  _  /mT7' 

+  (D227  +  4D226  +  D225)—  -  eA6  — 

h4  VL  y 


a88 


+ 


( 

3d127 


+  4D  15 
66 — 
2 


1 

h2 


+ 


(2D228  +  4D227 


+ 


1 

d226)“7  "  eA7 
h4 
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a12  =  -(v-m  +  d120  +  40  l\( —  I  ~  -  2  (d221  +  d220)~ 

1  66- )\L  /  h2  h4 


a78 


+  D126  +  4D 


1 

h2 


1 

2(D227  +  d226)~ 
h4 


1 

a13  =  d221  “ 
h4 


1 

a68  =  d226  — 
h4 


The  D's  and  A's  must  be  assigned  values  according  to  the  nodal  point 
represented.  In  the  neighborhood  of  the  jump  in  stiffness  all  the  D's  and 
A's  should  be  assigned  values  based  on  the  segment  (of  width  h)  centered 
on  the  nodal  points  except  for  Dgg  which  should  be  assigned  a  value  based 
on  the  segment  between  nodal  points. 

The  stability  equations  (BIO)  are  a  homogeneous  set,  and  to  get  a  nontriv¬ 
ial  solution  for  the  buckling  strain  the  determinant  of  the  matrix  must  vanish. 
The  buckling  strain  is  the  lowest  value  of  the  applied  strain  e  that  satis¬ 
fies  this  equation  for  various  values  of  m.  Computations  were  based  on  N  =  8. 
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TABLE  I.-  MATERIAL  PROPERTIES 


(a)  Graphite/epoxy 


Modulus  in  fiber  direction,  GPa  .  131 

Modulus  normal  to  fibers,  GPa . 13.0 

Shear  modulus,  GPa . 6.41 

Major  Poisson's  ratio . 0.38 

Density,  Mg/nP . 1.52 

Thickness,  mm/ply  ......  .  0.14 


(b)  Aluminum  honeycomb3 


Density, 

Major  transverse 

Minor  transverse 

kg/nP 

shear  modulus. 

shear  modulus, 

GPa 

GPa 

24 

0.117 

0.062 

32 

.  159 

.083 

48 

.262 

.138 

64 

.393 

.193 

80 

.524 

.262 

96 

.641 

.324 

112 

.786 

.372 

128 

.938 

.427 

144 

1,08 

.469 

160 

1.21 

.483 

aReference  4. 
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TABLE  II.-  DESIGN  LOAD  AND  MASS  FOR  SIMPLY  SUPPORTED  PANELS  76.2  cm  LONG 


Panel 

Honeycomb- 

core 
density, 
kg/m  ^ 

Theoretical 

buckling 

load, 

MN/m 

Theoretical 

buckling 

strain 

Mass/Plan 

kg/m 

area, 

2 

Calculated 

Measured 

Moderately  loaded 

48 

1.15 

0.0044 

6.93 

7.31 

160 

1.36 

.0051 

7.30 

8.39 

Lightly  loaded 

48 

.53 

.0046 

4.29 

5.73 

TABLE  III.-  BUCKLING  STRAINS  AND  MODES 


Panel 

length, 

cm 

Honeycomb- 

core 
density, 
kg/m  ^ 

Strain, 

e 

Buckles 

along  panel  length 

Present 

theory 

NASTRAN 

Experiment 

Present 

theory 

NASTRAN 

Experiment 

Moderately  loaded  panel 

design 

17.9 

48 

0.0047 

0.0047 

0.0045 

1 

4 

3 

26.7 

48 

.0043 

.0050 

.0048 

1 

al 

al 

30.7 

48 

.0043 

.0052 

.0042 

1 

a2 

al 

76.2 

48 

.0044 

.0040 

,0041 

2 

a2 

a2 

76,2 

160 

.0051 

.  0051 

.0058 

1 

al 

al 

Lightly  loaded  panel  design 

40.6 

48 

0.0046 

0.0052 

2 

7 

40.6 

48 

.0046 

.0050 

2 

7 

aIgnores  buckle  at  panel  ends  resulting  from  clamped-end  biaxial 
stress  field. 
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iables  used  in  minimum-mass  design  procedure 


8.53 


(a)  Moderately  loaded  design. 


Figure  3.-  Cross-sectional  dimensions  of  panel  designs  considered. 
(All  dimensions  are  in  centimeters.) 
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Figure  4.-  Comparison  of  analytical  results  for  compressive  buckling  strain  of 
moderately  loaded  sandwich-blade  stiffened  composite  panel  with  panel  buckle 
leng th. 


ure  5.-  Comparison  of  analytical  results  for  compressive  buckling  strain  of 
lightly  loaded  sandwich-blade  stiffened  composite  panel  with  panel  buckle 
length . 


sity,  kg/m 


buckling  strains  for  moderately  loaded  sandwich-blade  stiffened  composite  panel 


pressive  buckling  strains  for  lightly  loaded  sandwich-blade 
ffened  composite  panel.  Core  density,  48  kg/m^. 
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